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L  INTRODUCTION 

The  best  state-of-the-art  methods  for  the  de¬ 
tection  of  internal  discontinuities,  like  cracks  in 
solid  materials,  use  ultrasonic  waves.  Typically, 
quartz  or  barium  titanate  piezoelectric  transducers 
are  used  to  generate  the  acoustic  waves  and  to 
detect  the  reflections  from  the  various  dis¬ 
continuities  in  the  system.  A  critical  aspect  of  this 
type  of  system  is  the  coupling  of  the  transducer  to 
the  material  under  test.  Usually,  the  material  is 
immersed  in  a  liquid  that  increases  the  coupling 
efficiency.  Reflections  or  echoes  are  then  received 
from  the  metal-liquid  interfaces  at  the  sample 
surface  and  from  the  air-metal  interfaces  of  any 
internal  cracks  or  flaws.  By  observing  the  round- 
trip  times  of  these  echoes,  one  can  distinguish  the 
various  interface  reflections.  This  technique  is 
currently  in  use  in  industry,  but  it  has  some  serious 
drawbacks.  The  sample  must  be  placed  in  a 
continuous-flow  fixture  that  greatly  increases  the 
complexity  of  the  scanning  techniques  required. 
This  process  is  very  time-consuming  and  is  there¬ 
fore  not  conducive  to  production-line  testing. 

A  method  for  overcoming  the  difficulties  of  the 
conventional  technique  for  flaw  detection  has  been 
proposed  by  H.  P.  Kalmus  of  Harry  Diamond 
Laboratories  (HDL),  and  a  preliminary  evalua¬ 
tion  of  the  proposed  method  has  been  made  by 
Giglio.1  Basically,  the  technique  proposed  in¬ 
volves  replacing  the  generating  and  receiving  pie¬ 
zoelectric  transducers  in  presently  used  systems 
with  a  pair  of  lasers.  One  laser,  generally  a  high- 
power  pulse  laser,  is  used  to  generate  an  acoustic 
wave  in  the  sample  under  test,  and  the  second  laser 
records  the  motion  of  the  surface  of  the  sample 
caused  by  the  acoustic  waves  generated  by  the  first 
laser.  A  schematic  of  the  proposed  system  is  illus¬ 
trated  in  figure  1 .  The  fact  that  no  physical  contact 
with  the  sample  under  test  is  necessary  is  a  major 
advantage  over  the  presently  used  technique.  An 
optical  scanning  system  could  be  built  that  would 
provide  rapid  testing  of  similar  items,  and  entire 
sample  lots  could  be  tested  as  they  actually  come 
off  the  assembly  line. 

1  D  A.  Giglio,  Preliminary  Evaluation  of  a  Proposed  Laser-Operated  Flaw 
Detection  Technique,  Harry  Diamond  Laboratories,  HDL-TM  73-32  (December 
1973). 


This  report  presents  the  results  of  a  more 
extensive  evaluation  of  the  Kalmus  flaw-detection 
proposal.  A  description  is  given  of  the  acoustic- 
wave  generating  system  and  the  investigation  of 
the  parameters  that  affect  the  pulse  shape  and 
amplitude  of  acoustic  waves.  Also,  some  theory  on 
the  generation  of  acoustic  waves  by  the  laser  is 
presented.  The  investigation  of  a  readout  or  re¬ 
cording  system  that  uses  a  Michelson  interfero¬ 
meter  is  also  described.  Some  preliminary  flaw- 
detection  data  on  a  polished  metallic  sample  are 
discussed,  and  finally,  data  on  the  detection  of 
flaws  in  artillery  shells  are  presented. 

2.  ELECTRO-OPTICAL  RECEIVER 
SYSTEM 

An  acoustic  wave  generated  in  a  sample  under 
test  by  a  high-power  laser  pulse  propagates  as  a 
longitudinal  pressure  wave.  When  this  pressure 
wave  strikes  a  boundary,  the  surface  of  the  sample 
is  displaced  as  the  wave  is  reflected  by  the  boundary. 
The  surface  motion  at  the  boundary  can  be  detected 
with  a  piezoelectric  transducer  that  generates  an 
electric  signal  proportional  to  the  mechanical  im¬ 
pulse.  The  motion  of  the  surface  can  also  be 
detected  optically  by  using  a  laser  beam  in  con¬ 
junction  with  an  interferometer  and  some  type 
of  signal  processing.  Most  of  the  work  described  in 
this  paper  involves  the  optical  detection  system, 
although  some  interesting  data  were  also  obtained 
using  BaTi03  transducers. 

The  optical  detection  system  chosen  for  this 
work  was  a  Michelson  interferometer.  A  sche¬ 
matic  of  the  system  is  shown  in  figure  2.  The  laser 
beam  enters  the  interferometer,  where  part  of  it  is 
reflected  towards  a  fixed  mirror,  M2,  by  a  half- 
silvered  mirror,  The  remaining  part  of  the 
beam  proceeds  to  the  sample  under  test  (the 
sample  acts  as  a  movable  mirror).  The  reflections 
from  the  sample  and  mirror  are  then  recombined  as 
shown  in  the  figure.  If  the  path  difference  between 
the  two  beams  is  not  zero,  then  there  is  a  phase 
difference  between  the  two  beams  at  the  plane  of 
the  detector.  In  general,  this  leads  to  an  inter¬ 
ference  pattern  at  the  detector  plane  if  the  two 
reflecting  surfaces  are  not  perpendicular.  Any 
motion  of  the  sample  under  test  will  modulate  this 
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Figure  1.  Proposed  laser  flaw-detection  system. 
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Figure  2.  Simplified  Michelson  interferometer  receiver  system. 
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phase  and  therefore  change  the  portion  of  the 
interference  pattern  viewed  by  the  detector.  The 
sample  motion  will  therefore  cause  a  change  in  the 
light  intensity  seen  by  the  detector  and  a  corre¬ 
sponding  change  in  the  electrical  output. 

The  feasibility  of  this  type  of  detection  system 
for  measuring  small  displacements  or  vibrations 
has  been  demonstrated  by  many  authors.1 4  To 
evaluate  the  intensity  pattern  produced  at  the 
detector  plane,  the  alignment  configuration  shown 
in  figure  3  is  used.  For  simplicity,  it  is  assumed  that 
the  misalignment  angle  0  between  the  test  and 
reference  beams  is  in  one  plane  only.  Using 
complex  notation,  we  have  for  the  test  beam 

E,  =  K,E0ei(wt  ~  k<y«  +  2[x0  +  Ax(t)]|)jz  (j) 


and  for  the  reference  beam 


KjEoe'^1  "  kh  +  x  sin©  +  (yi  +  y0)coso]J^  ? 


where  E0  is  the  electric  field  amplitude  of  the  laser 
beam,  Kj  and  K2  are  proportionality  constants 
which  account  for  mirror  losses,  k  is  the  wave- 
number  2  7 r/A,  and  Ax(t)  is  the  time-varying  dis¬ 
placement  of  the  sample  under  test.  The  detector 
current  is  proportional  to  the  intensity,  I,  which  is 
given  by 

I  cc  E*  E 
or 

I  “  (K?  +  K \)  E0* 2  + 

2KjK2Eo  cos  [kx  sin©  +  k(y04-  yj) 
(cos©  —  l)  +  k(2y,  —  2x0  —  2Ax(t))]  ;(3a) 


5  D.  A.  Giglia,  Preliminary  Evaluation  af  a  Proposed  Laser  Operated  Flaw 
Detection  Technique,  Harry  Diamond  Laboratories,  HDL-TM  73-32  (December 
1973). 

^  H.  A.  Deferrari  and  F.  A.  Andrews,  Laser  Interferometric  Technique  for 
Measuring  Small-Order  Vibration  Displacements,  J  Acoustic.  Sac.  Am.,  39 
(1966),  979-980. 

^  S.  Sizgaric  and  A.  A.  Gundjian,  An  Optical  Hamadyne  Technique  far 
Measurement  af  A  mplitude  and  Phase  af  Subangstram  Ultrasonic  Vibrations, 
Proceedings  af  IEEE  (July  1969),  1313 

4  C.  R.  Brown,  G.  R.  Brown,  and  D.  H.  Niblett,  Measurement  af  Small  Strain 
Amplitudes  in  Internal  Friction  Experiments  by  Means  af  a  Laser  Interferometer,  J. 

Physics,  E:  Scientific  Instruments,  5  (1972),  966. 


then 

I  -  (K2 4  +  K 1)  E02  + 

2K,K2E02  cos  [C(x)  —  2kAx(t)]  ,  (3b) 

where 

C(x)  =  k[x  sin©  +  (y0  +  y,)(cos©  —  1)  + 
(2y,  -  2x0)]  .  (4) 

Now  let  Ax(t)  =  dsinnt,  where  d  is  the  ampli¬ 
tude  of  the  displacement  of  the  surface  and  n  is  the 
frequency  of  vibration.  The  signal  current  is  pro¬ 
portional  to  the  time-varying  component,  so 

is  cos  [a  sinnt  —  C(  X)]  ,  (5) 

where  a  is  the  modulation  depth  .  By  using 

trigonometric  identities,  we  obtain 

is  cos  (a  sinnt)  cosC(x)  4- 

sin  (a  sinnt)  sinC(x)  .  (6) 

Then,  using  the  Bessel  function  expansions,5 
cos(a  sinnt)  = 

oo 

J0(a)  4*  2  S  J2n(a)cos2nnt  (7a) 

n  =  1 

sin(a  sinnt)  = 

2  2  J2n  +  i(a)sin(2n  +  l)nt  ,  (7b) 

n  =  0 
we  get 

cos(a  sinnt)  =  J0(a)  +  2J2(a)  cos2nt  + 

2J4(a)cos4nt  +  .  .  .  (8a) 

and 

sin(a  sinnt)  = 

5  G.N.  Watson,  A  Treatise  on  the  Theory  of  Bessel  Functions,  Cambridge  Press 
(1952). 
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2J,(a)sinm  +  2J3(a)sin3ftt  + 
2J5(a)sin  5iit  +  .  .  . 
Expanding  each  term,  using 


(8b) 


J3(a)  =  fg  - 


a 


+ 


a 


768  30720 


J4(a)  = 


384  7680 


+ 


368640 


.  •  (10) 


Jn(z) 


we  get 
Jo(a) 


£  (~l)m  (1/2  Z)n+2m 
m  =  0  m  !  (n  +  m)  ! 


(9) 


If  a  <  1,  all  but  the  first  term  in  each  expansion 
can  be  neglected,  although  with  considerable  error 
for  a  close  to  one,  then 


1 —  +  —  + 

1  4  ^  64  ^ 


cosC(x) 


cos2fit 


+ 


Ji(a)  =  -J 

a2 

J2(a)  =  Y 


T6  + 


96 


+ 


384 


3072 


+  .  .  . 

■  +  .  .  . 


cos4flt  +  .  .  . 


+ 


(11) 


sinC(x)  |a  sinflt  +  sin3flt  +  .  .  JJ 
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Figure  3.  Setup  of  interferometer  for  analysis. 
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If  a  narrowband  detector  centered  at  fl  is 
used,  the  signal  current  is  given  by  is  a  sinC(x) 
sinfit.  Also,  if  a  <<  1,  a  wideband  detector  can 
be  used  and  the  harmonics  are  negligible.  For 
instance,  if  a  =  0.1,  the  sin3fit  term  is  propor¬ 
tional  to  0.00002  and  the  cos  2 fit  term  is  propor¬ 
tional  to  0.0025.  Also,  for  quadrature  detection 
C(x)  is  set  as  close  to  90  deg  as  possible,  so  that 
the  cosC(x)  term  becomes  zero  at  the  detector. 
Even  at  other  positions,  with  the  modulation  depth 
small,  the  cosine  term  has  no  time-varying 
component. 

The  actual  signal  current  will  now  be  given  by 
Is  K  s  2K,K2Eo  J  isdA  (12a) 


From  these  expressions,  some  conclusions  can 
be  drawn  about  the  detector  size  and  position  for 
maximum  signal.  However,  first  consider  the  static 
fringe  pattern,  which  can  be  plotted  from  equation 
(3b)  with  Ax(t)  =  0.  The  fringe  pattern  exhibits  a 
maximum  when  C(x)  =  nir,  n  even,  so  that  the 
distance  from  one  peak  to  another  is  given  by 
finding  the  change  in  x,  Sx,  required  to  cause  a 
change  in  C(x)  of  2  77,  or 

Ajkfxsin©  +  (y0  +  y^cos©  —  1)  + 

(2yj  —  2x0)]j  =  A[k(xsin©  +  ip)]  =  2n . 

Since  x  is  the  only  variable,  the  separation  in  the 
maxima  are  determined  by 


or 

Is  «  2K,K!SE5L  sm!11 

r  sinC(x)dx  ,  (12b) 

where  a  detector  that  is  rectangular  with  height  L 
and  length  x2  —  x1  is  assumed,  and  S  is  the 
sensitivity  of  the  detector.  Evaluating  the  integral: 

r*i 

sinC(x)dx  =  sink[xsin©  + 

(yo  +  y,)(cos©  -1)  +  (2y,  -  2x0)  ]dx  = 


sink(xsin©  +  \jj)  dx 


where 

=  (yo  +  y,)(cos0  -  1)  +  (2y,  -  2x0)  . 
Then 


r 


sinC(x)dx  = 


ksin© 


sin  < -r  [(x!  +  x2)  sin©  +  2ip] 


sin  |  |  [(x2  -  Xj)  sin©]  J  .  (13) 


k5x  sin©  =  277  or  Sx  = 


277 

ksin© 


X 

sin© 


This  means  a  dark  or  light  band  is  A/2sin©  wide, 
which  is  one-half  the  distance  between  peaks. 
Furthermore,  for  the  case  ^  =  0,  it  is  apparent 
that  there  is  a  maximum  at  x  =  0,  a  minimum  at 
x  =  A/2sin0,  a  maximum  at  A/sin0,  and  so  on. 


If  we  now  examine  equation  1 3,  we  see  that  the 
maximum  signal  is  observed  when 

k(x2  -  x,)  m77 

- 2 -  sln®  =  ~2~  m  °dd 

and 


sin© 


k(xi  +  x2) 

2 


+  kip 


mn 


m  odd; 


or,  with  m  =  1 , 

x2  “  x,  = 
and,  with  ijj  =  0, 

(Xi  +  X2) 

2 


X 

2sin© 


X 

4sin© 


Now  x2  —  x,  =  the  detector  length,  and 
xi  +  x2 

- r -  is  the  position  of  the  center  of  detector. 


9 


For  maximum  signal,  £  =  ^  ,  giving  a  sig¬ 

nal  current  of 

Is  °C  16  K,K2SEo  A  |  sinflt.  (16) 

Therefore,  if  the  modulation  depth  is  small,  the 
observed  signal  directly  reproduces  the  vibrational 
motion  of  the  sample.  In  experiments  where  a 
polished  sample  was  used,  the  fringe  pattern  was 
positioned  directly  on  the  detector  and  the  signal 
was  observed  in  this  manner.  Figure  4  shows  the 
variation  of  the  intensity  at  the  detector  for  various 
path  differences  between  the  test  and  reference 
beam. 

In  some  experiments,  diffuse  or  unpolished 
surfaces  were  used  as  the  test  object,  and  a 
focusing  lens  was  required  to  collect  the  diffusely 
scattered  radiation  from  the  test  sample.  The 


<ieti2  > 


PATH  DIFFERENCE 

Figure  4.  Intensity  distribution  at  detector  plane  for  setup  of  figure  3  with  6  =  0. 


Maximum  signal  is  therefore  observed  when  the 
detector  is  one  static  fringe  width  wide  (one  dark  or 
one  bright  band),  located  halfway  between  a 
minima  and  maxima. 

At  the  optimum  position  we  get 

Is  oc  2KJK2SEQL  sinHt 

2 


klsin0  \ 

ksin0  °1A1  \  2  J 


sin 


(14) 


or 


Is  oc  8  KjK2SEgA 

/  k£sin0 

TTd  .  Sm  V  2 

— r —  sinHt 


k£sin0 


•  (15) 
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roughness  of  the  test  sample  reduces  the  temporal 
and  spatial  coherence  of  the  reflected  radiation. 
This  effectively  reduces  the  signal-to- noise  ratio  at 
the  detector,  since  the  incoherent  part  of  the 
reflected  signal  acts  as  background  noise.  Also  in 
this  case,  the  spherical  aberration  introduced  by 
the  lens  created  a  circular  fringe  pattern  at  the 
detector.  Since  this  was  somewhat  awkward  to 
use,  the  misalignment  angle  was  adjusted  to  give  a 
parallel  fringe  pattern.  We  then  placed  a  variable 
slit  and  lens  in  front  of  the  detector  and  focused  the 
fringe  pattern  after  passing  through  the  slit  onto  the 
detector.  Using  this  technique  we  were  able  to  get 
high  sensitivities,  even  with  the  diffuse  surfaces. 


3.  LASER  SYSTEM  FOR  GENER¬ 
ATING  ACOUSTIC  WAVES 

3.1  Laser  Description 

A  commercial  ruby  laser  was  used  to 
generate  the  acoustic  waves.  The  laser  has  a  12.7 
cm  long  by  0.95  cm  diameter  Brewster/Brewster 
rod.  The  laser  was  run  in  the  Q-switched  mode 
with  a  Pockel’s  cell  Q-switch.  The  pulse  energy 
varied  between  0.9  and  1.4  J,  and  the  pulse  width 
was  50  ns  at  the  half-power  points.  The  power  level 
varied  from  1 8  to  28  MW.  The  unfocused  spot  size 
was  0.36  cm2,  which  gives  a  power  density  from  50 
to  78  MW/cm2.  Focused  power  densities  up  to 
1  GW/cm2  are  easily  obtained,  but  most  of  our 
experiments  involved  power  densities  under 
500  MW/cm2. 

3.2  Radiation  Pressure 

The  radiation  pressure  due  to  the  reflec¬ 
tion  of  the  optical  radiation  from  the  sample  for 
normal  incidence  can  easily  be  calculated.  The 
incident  intensity  at  the  sample  is 

I  =  Nhv,  (17) 

where  N  is  the  number  of  photons  per  second  per 
unit  area  incident  on  the  sample,  h  is  Planck’s 
constant,  and  v  is  the  frequency  of  the  radiation. 
For  a  reflection  coefficient,  R,  the  momentum 


transferred  to  the  sample  per  photon  upon  reflec¬ 
tion  for  a  plane  wave  is 

AM  =  (1  +  R)M  =  (1  +  R)  h/X  ,  (18) 

where  X  is  the  wavelength.  The  total  radiation 
pressure,  which  is  the  total  change  in  momentum 
per  unit  area  per  unit  time,  is 

P  =  (1  +  R)Nh/\  ,  (19) 


which  with  substitution  from  (17)  gives 
P  =  (1  +  R)h/X^=  (1  +  R)I/c,  (20) 

where  c  is  the  velocity  of  light. 

The  quantity  I/c  is  simply  the  energy  per 
unit  volume  of  the  laser  beam  ,  or  the  power 
density  divided  by  c.  For  a  power  density  of 

500  MW/cm2,  typical  for  our  experiments,  and 
a  reflection  coefficient  of  0.3,  appropriate  for  stain¬ 
less  steel  at  0.7  /xm,6 7 8 * * 11  the  radiation  pressure  is 
2.0  X  105  dyn/cm2. 

3.3  Thermal  Effects 

The  problem  of  thermoelastic  stress  pro¬ 
duction  in  solids  has  been  discussed  by  many 
authors.7-11  The  problem  consists  of  solving  two 
differential  equations.  First,  one  solves  the  heat- 
flow  equation  to  determine  the  temperature  profile 
in  the  material  generated  by  the  incoming  laser 
radiation.  Then  one  must  solve  the  equation  re¬ 
lating  the  thermoelastic  stress  to  the  temperature 
profile. 


6  M  Bam  and  F.  Wa\f,  Principles  of  Optics,  New  York,  Perga  mo n  Press 
(1959),  619 

7  J.  F.  Ready,  Effects  Due  to  Absorption  af  Laser  Radiation,  J.  Appl.  Phys., 
36,  2  (Februarv  1965),  462-468. 

8  J.  C.  Bushnell  and  D.  J.  McClaskey,  Thermae lastic  Stress  Production  in 
Solids,  J.  Appl.  Phys.,  39  (1968),  5541-5546. 

®  R.  M.  White,  Generation  of  Elastic  Waves  by  Transient  Surface  Heating,  J. 
Appl  Phys.,  34,  (1963),  3559-3567. 

*0  L.  S.  Gaumay,  Conversion  af  Electromagnetic  ta  Acoustic  Energy  by 
Surface  Heating,  J.  Acoustic.  Soc.  Am.,  40,  6  (1966),  1322-1330. 

11  T.  A.  Zaker,  Stress  Waves  Generated  by  Heat  Addition  in  an  Elastic  Solid,  J. 
Appl.  Mech.  (March  1965),  143-150. 
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Ready7  has  calculated  the  temperature 
distribution  in  a  semi-infinite  plane  slab,  assuming 
one-dimensional  heat  flow,  no  reradiation  of  en¬ 
ergy,  temperature-independent  thermal  properties, 
and  that  the  liquid  phase  can  be  neglected.  The 
heat  production  then  takes  the  form 

A(x,t)  =  I(t)  ae~x/S  ,  (21) 

where  8  is  the  optical  penetration  depth,  a  is  the 
absorption  coefficient,  and  I(t)  is  the  intensity  dis¬ 
tribution  of  the  laser  pulse  in  time.  If  one  assumes 
that 


equations  governing  the  propagation  of  stress 
waves  generated  in  an  isotropic  medium  due  to 
illumination  of  a  free  surface  by  a  laser  pulse  of 
energy  density  Q  were 


2(?,r)  = 

1/2 

~dr  ^y,r  ~~  ^  +  y^dy  + 

i 

1/2  ^ 

-fa  [0(y,r  +  £  -  y)]dy, 

*  <  s 

(23a) 

I(t)  =  Io  for  t  >  0  and  I(t)  =  0  for  t  <  0, 

one  obtains  the  solution  originally  derived  by 
Carslaw  and  Jaeger.12 


T(x,t)  =  (2I0a5/K)  (/ct)1/2ierfc  [x/2(/ct)1/2] 

-  (I0aS2e-x/5/K) 

+  (l0a52/2K)eKt/52_x/5erfc  j(Kt)1/2/5  - 


[x/2(Kt)1/2] 


(22) 


+  (l0a52/2K)eKt/52  +  x/S  erfc  ((«t)1/2/5 
[x/2(/ct)1/2]j  , 


+ 


where  K  is  the  thermal  conductivity  and  k  the 
thermal  diffusivity  K/pC,  with  p  the  density,  and  C 
the  specific  heat  of  the  material. 


Ready  then  applies  Duhamel’s  theorem12 
to  get  the  temperature  profile  for  a  laser  pulse  with 
arbitrary  shape.  For  appropriate  I(t)’s  repre¬ 
senting  real  laser  pulses,  the  temperature  must  be 
calculated  by  a  numerical  integration.  However,  to 
solve  the  thermoelastic  stress  equation,  we  need  a 
closed-form  solution  for  T(x,t). 


An  alternative  solution  has  been  given  by 
Bushnell  and  McCloskey.8 * * *  They  showed  that  the 

7  J  F.  Ready,  Effects  Due  to  Absorption  of  Laser  Radiation,  J.  Appl.  Phys., 
36,  2  (February  1965),  462-468. 

8  J.  C.  Bushnell  and  D.  J  McCloskey,  Thermoelastic  Stress  Production  in 

Solids,  J.  Appl.  Phys.,  39  (1968),  5541-5546. 

*2  H.  S.  Carslaw  andj.  C.  Jaeger,  Conduction  of  Heat  in  Solids,  Clarendon 

Press,  Oxford,  England  (1959),  80;  30-31. 


n&)  = 

1/2  X  It  [0(y,T  “  *  +  y)^dy  + 

1/2  -fa  [0(y,r  +  £  ~  y)]dy  - 

1/2  X  ~k  t0(y,r  -  *  ~  y^dy  ’ 

r  >  (23b) 

where  E  —  8a/ yQ  =  normalized  stress, 

£  =  x/5  =  normalized  distance, 

t  =  c0t /8  =  normalized  time, 

9  =  PC8( T  -  T0)/Q 

=  normalized  temperature, 

and  a,  8,  p,  c0,  C,  and  y  are  the  stress,  penetration 
depth,  density,  longitudinal  sound  speed,  specific 
heat,  and  Griineisen  constant  of  the  absorbing 
medium,  respectively,  and  y  is  a  dummy  integra¬ 
tion  variable. 

A  reasonable  description  of  0(£,r)  can  be 
obtained  by  assuming  that  the  absorbed  energy 
decays  exponentially  below  the  surface,  and  the 
pulse  width  is  short  enough  so  that  heat  conduc¬ 
tion  and  conversion  of  thermal  energy  to  acoustic 
energy  can  be  neglected.  Then 

Xr/r0 

g(y)  dy  ,  (24) 


12 


+4 


2 


where  g(y)  describes  the  laser  pulse  shape  in 
normalized  time  and  r0  =  c0t0/5,  the  ratio  of  the 
acoustic  penetration  depth  during  the  risetime  t0  of 
the  laser  pulse  to  the  optical  penetration  depth. 
Using  the  practical  envelope  description, 

g(y)  =  2(e-y  —  e_2y)  , 

we  get 

0(£,t)  =  e-i;  [l  —  2  exp  (— r/r0)  + 

exp  (— 2r/r0)]  .  (25) 

If  we  then  substitute  using  this  #(£,r)  in  the 
equations  for  E(£,r)  which  were  calculated  by 
Bushnell  and  McCloskey  in  a  manner  similar  to 
that  done  by  Giglio,1  we  get 

2(&r)  =  ~4  e~^+2r/ro>  + 


To  -  4 


ro  —  1 


e-(£+r/r0)  + 


(26a) 


1 


(To -2)  (To -1) 

1 

(To  +  1)  (t0  +  2) 


e_(r  +  4)  + 


r<  «  , 


and 


2(£r)  =  -j^-[e_2(r'^-e^  +  2r/ro>]  + 


1 


(To -2)  (To-1) 

2 


l_[V(r+o-e-(r-<n]  + 


=rt« 


T  >€ 


U+r)/r0  r)/r0 


]  , 

(26b) 


For  distances  large  compared  with  the  penetration 
depth,  and  making  the  substitution  7}  =  r  —  £, 
we  have 


*  Z).  Giglio,  Preliminary  Evaluation  of  a  Proposed  Laser-Operated  Flaw 
Detection  Technique,  Harry  Diamond  Laboratories,  HDL-TM-73-32  (December 
1973). 


Z(r,)  = 


tS-4 


p  2  - 


r0  —  1 


*~V/T0  — 


(t.-2)(r„-l)  ■  ">»  ■  <27a) 


and 


(r0+l)(r0  +  2)  ’  11  K  °  '  (27b) 

We  see  that  the  stress  wave  has  assumed  a 
traveling  form,  since 

>  _  (cot  -  x) 

V  =  T-  ? - g 


^  (wst  “  M)  ,  (28) 


where  o)s  and  ks  are  a  normalized  frequency  and 
wavenumber  dependent  on  the  penetration  depth 
S.  To  further  reduce  these  equations  we  need  to 
estimate  S. 

For  nonmagnetic  steel,  the  penetration 
depth  for  a  frequency  of  1014  Hz  is  —  10-6  cm. 
This  means,  with  a  50-ns  risetime  pulse,  that 
coto 

T°  ~  8 


=  2.9  X  104  . 


We  can  then  write 


T{n)  =  +  — ^ —  e_"/ro  -  4-  e~"/r°  - 


To 


X(V)  =  — 

r0 


,  n  >  0  ,  (29a) 

,  t]  <  0  (29b) 


The  above  equations  indicate  that  the  largest  value 
ofZ(rj)  is  given  by 

^max  =  l/ro  * 

or  in  terms  of  the  stress, 
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G’max 


=  yQ  =  yQS 

St2o  (c0to)2 


(30) 


where  Q  is  the  energy  density  of  the  illuminating 
pulse.  More  correctly,  we  must  replace  Q  by 
Q  (1  —  R),  which  corresponds  to  the  absorbed 
energy.  Thus,  the  peak  stress  is  given  by 


<W  =  Q  (1  -  R)  7~Ti  ’  (31) 

(c0t0)2 

where  we  have  taken  y  =  1  for  steel.  Also, 
Q  =  Poto/A,  where  P0  is  the  laser  peak  power  and 
A  the  beam  area.  Therefore, 


^max 


and  with 


Po 

A 


Po(l  ~R)8 

A  Cot0 


500 


MW 


cm 


R  =  0.3 

c0  =  6  X  103  m/s 
t0  =  50  X  10~9  s 
S  =  10  6  cm 


(32) 


we  get  (7max  =  2  X  105  dyn/cm2.  This  value  is 
about  twice  the  value  obtained  for  the  radiation 
pressure.  For  most  reasonable  values  of  R,  the 
thermal  pressure  is  somewhat  larger  than  the 
radiation  pressure. 


3.4  Blow-Off  Material  Pressure 


The  vaporization  of  material  from  the 
target  during  the  laser  pulse  can  contribute  strongly 
to  the  pressure  wave  generated  in  the  sample  under 
test.  The  material  removed  from  the  sample  during 
the  laser  pulse  imparts  a  momentum  to  the  sample 
that  is  equal  in  magnitude  and  opposite  in  sign  to 
that  of  the  blow-off  material.  This,  in  turn,  pro¬ 
duces  a  pressure  wave  in  the  sample.  The  depth  of 
material  removed  from  a  target  due  to  vaporiza¬ 
tion  by  Q-switched  laser  pulses  has  been  reported 


in  the  literature. 1 3  A  value  of  1 . 1  p,m  for  a  stainless- 
steel  target  and  a  power  density  of  109  W/cm2  is 
typical.  In  our  experiments,  a  power  density  of 
5  X  108  W/cm2  produces  the  largest  signal;  the 
depth  of  material  removed  per  shot  is  there¬ 
fore  —  0.5  jam. 


A  simple  estimate  of  the  pressure  exerted 
on  the  sample  by  the  vaporized  material  can  be 
made  if  the  volume  of  material  removed  and  the 
velocity  of  this  material  is  known.  Weichel  and 
Avizonis14  have  measured  the  plasma  plume  veloci¬ 
ty  observed  when  illuminating  a  carbon  target.  Their 
measurements  showed  that  the  initial  velocity  of 
the  plume  was  approximately  the  same  for  any 
energy  density  which  produced  vaporization.  At 
higher  energy  densities,  the  plume  velocity  in¬ 
creases  from  the  initial  value  due  to  the  absorption 
of  the  laser  radiation  and  subsequent  heating  of  the 
plasma.  The  energy  densities  used  in  our  experi¬ 
ments  were  only  slightly  above  the  point  where 
vaporization  occurs,  so  the  initial  vaporization 
velocity  will  be  used  in  our  calculations. 


Weichel  and  Avizonis  found  that  the 
initial  velocity  approached  the  normal  vaporiza¬ 
tion  velocity,  which  can  be  estimated  from 

Vi  -  (f  CTS^  ,  (33) 

where  C  is  the  specific  heat  and  Ts  is  the  sub¬ 
limation  temperature  of  the  target.  For  steel,  the 
sublimation  temperature  will  be  taken  to  be 
'v,  3000°K,  so  that,  with  C  =  0.107  cal/g°K,  we 
get  Vi  =  2  X  105  cm/s.  The  volume  of  material 
removed  at  500  MW/  cm2  can  be  calculated  from 
the  vaporized  depth  (0.5  fxm)  and  the  spot  area 
(0.055  cm2).  The  momentum  of  the  blow-off 
material  is  then  given  by 

Mb  =  V  V-P  ,  (34) 

where  V  is  the  volume  and  p  the  target  density. 

^  J.  F.  Ready,  Effects  of  High- Power  Laser  Radiation,  New  York,  Academic 
Press  (1971), 111,  161-206. 

1 4  H  Weich  el  and  P.  V.  A  vizonis.  Expansion  Rates  of  the  L  uminous  Front  of a 
Laser-Produced  Plasma,  Appl.  Phys.  Lett.,  9,  9  (1966),  334-337. 


14 


With  V  =  2.75  X  1CT6  cm3 

Vj  =  2  X  105  cm/s  ,  and 
p  =  7.75  g/cm3  , 

Mb  would  be  4  dyn-s.  To  obtain  the  pressure,  we 
divide  by  the  pulse  width  and  area.  The  blow-off 
pressure  is  then  1.5  X  109  dyn/cm2.  Measure¬ 
ments  of  the  momentum  transferred  to  a  sample 
from  surface  vaporization  have  been  reported  in 
the  literature.  Gregg  and  Thomas15  reported  a  mo¬ 
mentum  transfer  of  5/dyn-s/J  for  an  aluminum 
target  using  a  7.5-ns  ruby  laser  pulse  with  peak 
power  density  in  the  range  of  500  MW/cm2.  The 
pulse  energy  was  0. 1 1  J,  giving  a  total  impulse  to  the 
target  of  0.5  5  dyn-s.  Neuman,16  in  a  similar  experi¬ 
ment,  measured  impulses  of  0.17  dyn-s  for  both 
aluminum  and  cold  rolled  steel  targets.  The  peak 
power  density  was  800  MW/cm2,  the  pulse  width 
was  50  ns,  and  the  total  energy  was  0.3  J. 

To  convert  these  impulse  values  to  pres¬ 
sure  at  the  surface  of  the  sample,  we  divide  by  the 
pulse  length  and  spot  size.  For  comparison  to  the 
values  calculated  previously  for  radiation  and  ther¬ 
mal  pressure  we  also  normalize  to  the  same  total 
energy.  With  a  1-J,  50-ns  pulse  and  a  spot  size  of 
0.4  cm2,  the  pressures  for  the  cases  measured  by 
Gregg  and  Thomas  and  Neuman  would  be 
2.5  X  109  dyn/cm2  and  2.8  X  108  dyn/cm2, 
respectively.  The  values  are  much  higher  than  the 
estimated  radiation  pressures  and  in  reasonable 
agreement  with  the  value  obtained  from  the  esti¬ 
mate  based  on  material  removal  alone. 

A  significant  difference  exists  between 
our  experiments  and  those  done  by  Gregg  and 
Thomas  and  Neuman.  Their  measurements  were 
performed  in  a  vacuum  whereas  ours  were  in  air. 
The  presence  of  air  at  the  focal  position  of  the  laser 
pulse  could  increase  the  degree  of  plasma  forma¬ 
tion  due  to  the  higher  density  of  atoms  in  the 
plasma  volume.  These  atoms  can  contribute  to  the 
cascade  process  involving  free  electrons  that  takes 
place  during  plasma  formation.  Once  a  significant 


I5  D  W  Gregg  and  S.  J.  Thomas,  Momentum  Transfer  Produced  by  Focused 
Laser  Giant  Pulses,  J.  Appl.  Phys.,  37  (1966),  2787. 

F.  Neuman,  Momentum  Transfer  and  Cratering  Effects  Produced  by  Giant 
Laser  Pulses,  Appl  Phys.  Lett.,  4  (1964),  167. 


amount  of  material  is  vaporized,  the  plasma  will 
shield  the  sample  from  the  later  parts  of  the  laser 
pulse.  The  absorption  of  the  laser  pulse  by  the 
plasma  will  increase  the  plasma  density  and  reduce 
the  measured  momentum  transfer.  Therefore,  we 
might  expect  slightly  lower  numbers  for  our  experi¬ 
ments.  As  indicated  by  Gregg  and  Thomas15  and  in 
Neuman,16  the  effect  of  plasma  formation  becomes 
important  at  power  levels  greater  than  1  GW/cm2, 
even  in  a  vacuum.  As  the  power  density  increases 
beyond  this  point  the  impulse  delivered  to  the 
sample  decreases.  The  power  level  of  peak  mo¬ 
mentum  transfer  is  probably  lower  in  a  non- 
evacuated  system. 

If  the  blow-off  material  pressure  is  com¬ 
pared  to  the  thermal  and  radiation  pressures  as 
calculated  at  500  MW/cm2,  the  blow-off  pressure 
is  between  three  and  four  orders  of  magnitude 
higher  than  the  thermal  or  radiation  pressure.  The 
blow-off  material  pressure  is  easily  the  dominant 
pressure  at  500  MW/cm2.  However,  if  the  power 
density  is  reduced  by  10,  the  thermal  and  radiation 
pressures  drop  linearly  while  the  blow-off  pres¬ 
sure,  which  depends  on  the  amount  of  ionized 
material  ejected  from  the  surface,  drops  much 
faster.  The  experimental  observation  of  this  effect 
will  be  discussed  in  a  later  section  in  some  detail. 


3.5  Estimation  of  Expected  Surface 
Motion 

To  calculate  the  motion  of  the  surface  of 
the  sample  caused  by  the  impulsive  force  delivered 
by  the  laser  pulse,  a  few  assumptions  are  necessary. 
First,  we  assume  that  the  pressure  is  uniformly 
exerted  on  the  sample  over  the  half-power  pulse 
width  of  the  laser  pulse.  Secondly,  we  assume  that 
the  strained  material  depth  is  equal  to  the  depth  of 
penetration  of  the  resulting  acoustic  wave  during 
the  laser  pulse  duration.  We  also  assume  that  the 
acoustic  wave  generated  in  this  manner  is  a  plane 
wave.  This  implies  that  the  stress  and  strain  in  the 
sample  are  related  by 

exx  =  P  ,  (35) 

'-'ii 
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where  exx  =  strain, 

P  =  pressure,  and 

Cn  =  compliance  coefficient. 

Let  us  assume  that  P  =  108  dyn/cm2, 
due  to  blow-off  for  a  power  density  of 500  MW/cm2. 
This  is  a  conservative  estimate,  since  the  value  is 
somewhat  lower  than  that  estimated  from  the 
literature  discussed  in  the  previous  section.  This 
gives,  with  Cn  =  2.4  X  1012  dyn/cm2  for 
steel,17  a  value  for  exx  of  4.2  X  1CT5. 

To  calculate  the  motion  of  the  front  sur¬ 
face  of  the  sample,  we  multiply  the  strain  by  the 
penetration  depth  of  the  stress  pulse,  which  is  just 

Pd  =  c0t  ,  (36) 

where  c0  =  5.92  X  105  cm/s,  the  longitudinal 
sound  speed  for  steel,6  and  t  =  50  X  10~9  s,so 
that 

Pd  =  3.0  X  10-2  cm 

Therefore, 

A^  =  Pd  exx  =  1.2  X  10-6  cm  ,  or 

A^  =  0.012  jam  for  F  =  108  dyn/cm2 

Assuming  the  minimum  blow-off  material 
pressure  from  the  literature  of  F  =  2.8  X 
108  dyn/cm2,  we  would  get  a  Af  =  0.034  ^um.  It 
must  be  stressed  that  this  motion  is  for  the  front 
surface  only.  It  also  implies  that  all  momentum 
transferred  to  the  sample  is  contained  in  a  narrow 
acoustic  pulse.  Divergence  or  beam  spread  of  the 
acoustic  wave  and  attenuation  in  the  material 
results  in  a  smaller  A  2  at  the  point  of  observation. 
In  our  readout  system,  we  used  both  a  NeHe  laser 
and  an  Argon  laser.  The  operating  wavelengths 
were  632.8  and  514.5  nm,  respectively.  With 
these  lasers  and  a  maximum  surface  motion  of 
0.034  jam,  the  maximum  signal  for  the  direct 
wave  would  be  in  the  A/1 5  range.  The  intensity  of 


Born  and  F.  Wolf,  Principles  of  Optics,  New  York,  Pergamon  Press 
(1959),  619. 

”S.  P.  Clark,  ed.,  Handbook  of  Physical  Constants,  Geological  Soc.  of 
America  (1966). 


the  acoustic  wave  scattered  from  a  flaw  is  less,  of 
course,  since  only  part  of  the  acoustic  wave  is 
scattered.  The  observed  signal  from  a  flaw  is  a 
complicated  function  of  the  flaw  size,  position,  and 
geometry.  A  detailed  theoretical  evaluation  of 
these  scattering  effects  is  beyond  the  scope  of  this 
report. 

3.6  Analysis  of  Plasma  Effects 

A  high-density  plasma  is  created  when 
the  laser  beam  vaporizes  the  surface  of  the  test 
sample.  The  density  and  velocity  of  this  plasma 
can  significantly  affect  the  amount  of  radiation 
reaching  the  surface  and  the  total  amount  of 
material  removed.  Ready13  presents  a  comprehen¬ 
sive  review  of  the  experimental  and  theoretical 
work  which  has  been  done  in  this  area.  Estimates 
of  the  plasma  density  and  velocity  for  our  experi¬ 
ments  are  based  on  his  book. 

The  plasma  electron  density  determines 
the  frequency  of  radiation  which  can  be  trans¬ 
mitted  through  the  plasma.  The  plasma  frequency 
isgivenbyvp  =  8.9  X  103  ne1/2,  where ne  is  the 
electron  density.  If  v  >  vp  the  incident  radiation 
is  reflected  by  the  plasma,  and  if v  <  vp  the  radia¬ 
tion  is  transmitted.  For  the  ruby  wavelength,  the 
radiation  would  be  reflected  if  ne  >  2.4  X 
1021  cm-3.  The  electron  density  can  approach 
1019  cm-3  for  short  times  with  lasers  similar  to 
ours.  The  very  high  velocity  of  plasma  expansion, 
typically  1 06  to  1 07  cm/s,  would  reduce  the  electron 
density  quickly.  It  is  obvious  that  the  electron 
density  does  not  approach  that  needed  to  reflect 
the  ruby  light  in  our  case. 

Experimental  evidence13  shows  that  the 
plasma  is  accelerated  during  the  laser  pulse.  This 
means  that  although  the  vapor  leaves  the  surface  at 
a  velocity  commensurate  with  the  normal  vapori¬ 
zation  velocity  due  to  boiling,  the  velocity  can 
increase  by  orders  of  magnitude  after  leaving  the 
surface.  This  increase  in  velocity  is  due  to  the 
absorption  of  photons  by  free  electrons  in  the 
inverse  Bremsstrahlung  process.  The  electrons  are 
raised  to  a  higher  state  in  the  continuum.  Weichel 

^  J.  F.  Ready,  Effects  of  High-Power  Laser  Radiation,  New  York,  Academic 
Press  (1971),  161-206,  111. 
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and  Avizonis14  have  shown  that  for  an  energy 
density  of  70  J/cm2,  the  plasma  expansion  velocity 
can  increase  by  as  much  as  five  times.  This 
absorption  would  reduce  the  amount  of  laser  light 
reaching  the  sample.  A  rough  estimate  of  the 
plasma  absorption  constant  can  be  made  by  using 
the  formula 

Kv  =  3.69  X  108  (ZV/T1/2v3) 

(l  -  e-hv/KBT)  ,  (37) 

where  Z  is  the  average  charge  per  ion,  n*  the 
density  of  ions,  T  the  temperature  of  the  plasma, 
and  Kb  the  Boltzmann  constant.  Measured  plasma 
temperatures  for  a  power  density  of  5  X 
108  W/cm2  are  about  5  eV.  We  will  assume  that 
ne  is  1019  cm-3.  For  iron,  which  has  26  electrons, 
the  outside  shell  electrons  are  very  loosely  bound. 
Spectra  of  Fe(XII)  and  Fe(XIII)  have  been  ob¬ 
served  in  a  plasma  generated  by  a  laser  pulse 
striking  an  iron  target.  Experiments  also  indicate 
that  the  highly  ionized  states  appear  early  in  the 
plasma,  whereas  the  lower  ionized  states  appear 
later.  Therefore,  we  will  assume  that  Z  is  10.  This 
would  make  nj  =  1018  cm-3  for  ne  =  1019cm-3. 
Kv  then  equals  6.3  cm-1.  This  estimate  was  made 
from  data  that  were  taken  in  a  vacuum.13  Our 
measurements  were  made  at  normal  pressure;  the 
presence  of  air  at  the  focal  position  could  increase 
the  plasma  electron  density  because  of  the  extra 
ions  available  to  feed  the  cascade  process.  There¬ 
fore,  if  the  plasma  was  1  mm  thick,  about  50  per¬ 
cent  of  the  remaining  part  of  the  laser  pulse  would 
be  absorbed.  This  absorption  of  the  laser  light 
would  increase  the  plasma  temperature,  which 
would  in  turn  heat  the  sample  until  it  dissipated. 
Thus,  if  you  were  monitoring  the  energy  transfer  to 
the  sample,  you  would  first  see  an  energy  deposi¬ 
tion  which  followed  the  laser  pulse  until  vaporiza¬ 
tion  occurred.  At  this  point,  depending  on  the 
density  of  this  resulting  plasma,  there  would  be  a 
decrease  in  the  absorbed  energy  while  the  plasma 
was  heated.  The  plasma  would  then  transfer  heat 
to  the  sample  by  radiation  and  conduction  before  it 
completely  dissipated. 


13  J-  F.  Ready,  Effects  of  High-Power  Laser  Radiation,  New  York,  Academic 
Press  (1971),  111,  161-206. 

H  Weicheland  P.  V.  Avizonis,  Expansion  Rates  of  the  Luminous  Front  of  a 
Laser-Produced  Plasma,  Appl.  Phys.  Lett.,  9,  9  (1966),  334-337. 


4.  EXPERIMENTAL  SETUP  AND 
RESULTS 

4. 1  Transducer  Experiments 

Experiments  were  conducted  to  charac¬ 
terize  the  acoustic  waves  generated  in  a  steel 
sample  by  a  high-powered  laser  beam.  A  sche¬ 
matic  of  the  experimental  setup  is  shown  in  figure 

5 .  The  sample  used  most  often  was  a  steel  cylinder 
5  cm  in  diameter  and  2.5  cm  thick.  A  BaTi03 
transducer  was  placed  on  the  back  face  of  the 
sample  under  test.  The  transducer  had  a  diameter 
of  6.35  mm  and  a  thickness  of  3.17  mm,  with  a 
thickness-mode  resonance  of  800  kHz  and  a 
radial-mode  resonance  of  400  kHz.  The  trans¬ 
ducer  was  placed  against  the  back  of  the  steel 
sample  with  silicone  grease,  and  pressure  was 
applied  to  its  back  face  by  a  rubber-covered  screw. 
The  sample  was  grounded  to  the  outer  shield  of  a 
50-ft  cable  while  the  center  conductor  of  the  cable 
was  conected  to  the  back  face  of  the  transducer. 
The  electrical  load  was  supplied  by  a  termination 
at  the  oscilloscope. 

The  loading  conditions  on  the  transducer 
have  a  large  effect  on  the  Q  of  the  transducer  and 
therefore  its  bandwidth.  Huester  and  Bolt18  have 
shown  that  by  loading  a  transducer  with  steel  and 
using  a  rubber  backing,  the  mechanical  Q  of  the 
transducer  can  be  reduced  drastically,  and  by 
terminating  the  transducer  in  a  low  impedance,  the 
electrical  Q  can  also  be  reduced.  Figure  6  shows 
the  effect  of  changing  the  load  impedance  from 
1  Mn  to  50  ft.  In  figure  6(a)  the  signal  received 
by  the  transducer  terminated  in  1  Mft  is  shown  for  a 
laser  pulse  of  0.9  J  or  50  MW/cm2.  The  laser  pulse 
is  incident  on  the  front  face  of  the  sample  near  the 
beginning  of  the  trace.  Its  position  is  indicated  by 
the  ringing  due  to  pickup  in  the  input  cable  caused 
by  the  high  voltage  discharge  of  the  laser.  The 
delay  in  the  signal  from  the  laser  pulse  to  the 
acoustic  pulse  is  due  to  the  transit  time  of  the 
acoustic  wave  through  the  sample.  In  figure  6(b)  the 
same  pulse  is  shown  with  the  cable  terminated  with 
a  50-ft  load.  Notice  that  the  signal  level  is 
decreased  and  that  the  pulses  have  become  much 

18  T.  F.  Huester  and  R.  H.  Bolt,  Sonics,  New  York,  John  Wiley  and  Sons  Inc. 
(1955),  106-109. 
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Figure  5.  Experimental  setup  for  transducer  measurements. 
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Figure  6.  Received  signals  using  BaTi03  transducer  for  different  loading  conditions:  (a)  1-M  ft  termina¬ 
tion,  (b)  and  (c)  50-  a  termination,  (d)  incident  laser  pulse. 
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narrower.  Figure  6(c)  shows  the  first  two  pulses  in 
figure  6(b)  on  an  expanded  time  s  cale  ( 1 00  ns/div). 
The  first  pulse  is  a  good  reproduction  of  the  inci¬ 
dent  laser  pulse,  a  trace  of  which  is  shown  in  figure 
6(d).  The  second  pulse  is  the  reflection  of  the 
acoustic  wave  from  the  back  face  of  the  transducer. 
The  reason  for  the  increase  in  apparent  signal  level 
for  the  reflection  is  not  clear,  but  we  suspect  it  is 
caused  by  the  different  loading  conditions  at  the 
back  surface.  Since  the  pulse  shape  is  a  good 
reproduction  of  the  laser  pulse,  it  is  clear  that  the 
effective  bandwidth  of  the  transducer  has  been 
greatly  increased  by  loading  it  with  50  11,  and  it  is 
also  evident  that  the  ringing  in  the  transducer  seen 
in  the  later  part  of  figure  6(a)  has  been  almost 
completely  eliminated.  The  transducer  signal 
observed  using  a  50-H  load  indicates  that,  with  a 
50-MW/cm2  pulse  of  0.9  J  incident  on  a  piece  of 
steel,  the  resulting  acoustic  wave  shape  closely  ap¬ 
proximates  that  of  the  incident  laser  pulse.  This 
will  be  discussed  in  more  detail  later,  along  with 
some  other  acoustic  wave  characteristics. 

4.2  Direct  Transmission  Interferometer 
Experiments 

The  electro-optical  receiver  system  was 
discussed  in  some  detail  in  a  previous  section.  In 
this  section  the  experimental  results  for  the  de¬ 
tection  of  transmitted  acoustic  waves  using  the 
interferometer  will  be  discussed.  In  all  cases  the 


sample  was  a  steel  cylinder  2.5  cm  thick  and  5  cm 
in  diameter,  with  the  back  face  ground  and 
polished.  The  polished  side  was  used  as  one  mirror 
of  the  interferometer.  The  experimental  setup  is 
shown  in  figure  7.  A  narrowband  filter  precedes  the 
detector  amplifier  assembly  which  is  a  commercial 
silicon  avalanche  detector  module.  The  detector 
amplifier’s  maximum  sensitivity  is  3.7  X  104  V/W 
at  900  nm  and  slightly  less  at  630  nm.  The  out¬ 
put  from  this  assembly  was  split  into  two  parts 
as  shown  in  the  figure.  The  high-frequency  output 
was  fed  directly  to  a  500-MHz  bandwidth  oscillo¬ 
scope,  while  the  dc  output  was  fed  to  another 
oscilloscope  for  monitoring  the  fringe  position  at 
the  detector.  A  micrometer  adjustment  on  the 
interferometer  allowed  positioning  of  the  fringe 
pattern  so  that  the  detector  was  at  the  most  sensi¬ 
tive  position,  that  is,  halfway  between  a  minima 
and  maxima.  The  interferometer  was  typically  set 
up  so  that  the  difference  from  a  minima  to  a 
maxima  was  between  1  and  1.5  V.  A  complete 
swing  from  a  minima  to  a  maxima  corresponds  to  a 
phase  difference  of  180  deg  or  a  movement  at  the 
sample  of  1/4  of  a  wavelength,  or 0. 158  jam  when 
using  the  HeNe  laser.  Experiments  were  con¬ 
ducted  to  characterize  acoustic  wave  shape,  ampli¬ 
tude,  and  beam  spread  as  a  function  of  generating 
laser  parameters. 

A  typical  signal  from  the  electro-optical 
receiving  system  is  shown  in  figure  8.  The  initial 
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Figure  7.  Experimental  setup  of  direct  transmission  interferometer  measurements. 
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part  of  the  pulse  corresponds  to  the  acoustic  wave 
generated  by  the  expulsion  of  material  from  the 
surface,  while  the  slow  rising  wave  that  follows  is 
due  to  heating  effects  in  the  material.  To  measure 
the  angular  spread  of  the  acoustic  wave  in  the 
sample  the  interferometer  was  set  up  with  the 
generating  spot  and  readout  spot  centered  on 
opposite  sides  of  the  sample.  By  translating  the 
focusing  lens  shown  in  figure  7,  the  generating  spot 
could  be  moved  across  the  sample.  By  varying  the 
distance  between  the  focusing  lens  and  the  sample, 
the  generating  spot  size  could  be  varied.  Figure  9 
shows  a  plot  of  the  initial  peak  value  for  two  gener¬ 
ating  spot  sizes  as  a  function  of  position.  The  spot 
diameters  in  the  scan  direction  are  5.6  and  2.4  mm. 
The  widths  of  the  acoustic  beams  to  the  half-peak 
points  for  the  two  traces  are  8.1  and  13.8  mm, 
respectively.  This  is  what  one  would  expect,  since 
the  smaller  the  source,  the  closer  it  approximates  a 
point  source,  and  the  more  spherical  the  wave 
fronts  become.  This  indicates  that  the  smaller  the 
generating  spot,  the  larger  the  angular  spread  and 
the  faster  a  sample  could  be  scanned,  but  at  the 
expense  of  requiring  a  greater  detection  system 
sensitivity. 

Other  parameters  of  interest  are  the 
acoustic  wave  shape  and  amplitude.  With  the 
generating  spot  centered  with  respect  to  the  read¬ 
out  spot,  the  power  density  at  the  sample  was 
varied  in  two  ways.  In  one  case,  the  ruby  laser  was 
filtered  to  vary  the  energy  at  the  sample  while  the 
spot  size  was  held  constant.  This  was  done  for  two 
spot  sizes.  The  curves  in  figure  10  were  obtained 


with  the  0.3  cm2  unfocused  incident  beam.  The 
initial  peak  level  is  the  maximum  level  of  the  initial 
peak  shown  in  figure  8,  while  the  slow  peak  is  the 
maximum  level  of  the  slower  wave  shown  in  the 
same  figure.  The  fast  wave  increases  quite  drastic¬ 
ally  at  45  MW/cm2;  this  power  density  corres¬ 
ponds  to  the  point  where  material  is  starting  to  be 
removed  from  the  sample.  A  similar  series  of  data 
points  was  obtained  with  an  incident  beam  spot 
area  of  0.06  cm2.  The  maximum  power  density  in 
this  case  was  510  MW/cm2.  Figure  11  shows 
some  initial  pulse  signals  for  various  power 
densities.  Note  that  for  power  densities  above 
150  MW/cm2,  the  signal  has  a  different  shape. 
The  fast  rise  of  the  pulse  is  the  same,  but  the  signal 
does  not  return  to  the  baseline  level.  Figure  12 
shows  a  plot  of  the  initial  peak  signal  level  and  the 
signal  level  at  the  third  graticule  line  or  600-ns 
point  plotted  against  power  density  for  these  data. 

In  the  second  case  a  set  of  data  was  taken 
with  constant  energy  but  varying  spot  sizes.  Figure 
13  shows  some  typical  signal  returns  for  power 
densities  from  500  to  1000  MW/cm2.  Notice  that 
for  some  of  the  higher  levels,  the  initial  peak  is  not 
easily  visible  and  the  signal  level  actually  rises  to 
higher  values  later.  Figure  Mis  a  plot  of  these  data 
in  the  same  format  as  figure  12.  Itis  evident  that  the 
signal  level  for  both  the  initial  peak  and  600-ns 
point  drops  steadily  above  600  MW/cm2.  As  dis¬ 
cussed  previously,  this  is  due  to  the  formation  of  a 
plasma  at  these  high-power  densities;  the  plasma 
becomes  highly  absorbing  to  the  laser  radiation 
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Figure  8.  Typical  acoustic  waves  as  detected  with  interferometer. 
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Figure  9.  Acoustic  wave  beam  width  at  back  face  of  sample  for  two  generating  spot  diameters.  The  mea¬ 
sured  beam  widths  were  8.1  and  13.8  mm  for  5.6-  and  2-mm  diameter  spots,  respectively. 
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Figure  10.  Acoustic  wave  signal  level  versus  power  density  for  a  constant  generating  spot  size. 
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Figure  11.  Oscilloscope  photographs  of  received  acoustic  waves  for  various  power  densities,  keeping 
generating  spot  size  constant  The  vertical  scales  are  50  mV/ div  for  all  but  the  66  MW/cm2 
case,  when  it  is  20  mV/div.  Horizontal  scales  are  all  200  ns/div. 
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Figure  12.  Observed  acoustic  wave  amplitude  versus  power  density  for  two  different  positions  on  oscillo¬ 
scope  traces.  Data  were  obtained  using  reduced  spot  size  from  that  for  data  of  figure  10. 
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Figure  13.  Oscilloscope  photographs  of  acoustic  waves  obtained  by  keeping  pulse  energy  constant  and 
varying  generating  spot  size.  The  scales  are  20  mV/div  and  200  ns/div. 
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Figure  14.  Signal  amplitude  versus  power  density  for  two  positions  as  in  Figure  12. 
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and  thereby  shields  the  sample  from  the  latter  part 
of  the  laser  pulse.  The  plasma  is  heated  while  ab¬ 
sorbing  the  laser  pulse,  and  it  in  turn  transfers  some 
of  the  heat  to  the  sample  until  it  dissipates,  giving 
the  high  levels  at  the  600-ns  point. 

Figures  10,  12,  and  14  cannot  be  com¬ 
pared  directly  because  of  the  various  spot  sizes. 
For  a  more  detailed  evaluation  of  the  effects  of 
power  density  on  acoustic  wave  amplitude,  a  much 
higher  power  laser  would  be  needed  so  that  an 
entire  range  of  power  densities  of  interest  could  be 
scanned  without  varying  the  spot  size. 

A  comparison  between  calculated  and 
measured  surface  motion  is  difficult  to  make.  Data 
in  the  literature  indicate  that  the  maximum  impulse 
is  delivered  to  a  target  for  power  levels  of  about 
500  MW/cm2.14  These  data  were  taken  by 
varying  laser  energy  while  keeping  spot  size 
constant  at  0.03  cm2.  A  typical  surface  motion  for 
this  power  density  and  spot  size  in  our  experiments 
was  0.02  fim ,  compared  to  a  calculated  value  at 
the  front  surface  of  0.034  fim  (sect.  3.9).  By 
comparing  the  amplitude  of  the  directly  transmitted 
wave  with  that  of  the  wave  which  traverses  the 
sample  three  times,  the  one-way  transmission  was 
calculated  to  be  56  percent.  This  would  account 
for  the  difference  in  the  calculated  front  surface 
motion  and  the  observed  surface  motion  at  the 
back  face. 

The  maximum  value  of  the  initial  peak 
signal  observed  in  our  experiments  was  for  a  power 
density  of  70  MW/cm2.  This  is  a  lower  power 
density  than  reported  in  the  literature  but  with  the 
same  incident  energy  and  a  larger  spot  size.  It  is  the 
maximum  power  density  that  we  could  obtain  with 
this  spot  size.  Increasing  the  power  density  by  re¬ 
ducing  the  spot  size  reduced  the  signal  because  of  a 
combination  of  two  factors.  First,  the  smaller  spot 
size  creates  a  less  collimated  acoustic  beam  in  the 
sample.  Second,  the  higher  power  density  increases 
the  amount  of  plasma  formation,  which  means  that 
a  smaller  portion  of  the  energy  is  used  to  create  the 
acoustic  wave.  The  surface  displacement  at  the 
70  MW/cm2  power  density  was  0.03  am  at  the 

14  H.  Weichel  and  P.  V.  Avizonis.  Expansion  Rates  of the  Luminous  From  of a 
Laser  Produced  Plasma,  Appl.  Phvs.  Lett.,  9.  9,  (1966),  334-337. 


back  face.  Considering  the  assumptions  and  ap¬ 
proximations  made  in  calculating  the  surface  dis¬ 
placement,  along  with  the  range  of  data  in  the 
literature,  the  values  obtained  are  in  good  agree¬ 
ment  with  calculated  values. 

4.3  Reflection  and  Flaw-Detection 

Experiments 

In  order  to  demonstrate  the  feasibility  of  a 
remote  flaw-detection  system  where  the  acoustic 
waves  are  both  generated  and  detected  by  laser 
beams,  the  experimental  setup  shown  in  figure  15 
was  used.  To  increase  the  dynamic  range  and 
sensitivity  of  the  readout  system  over  the  system 
shown  in  figure  7,  an  Argon  laser  operating  with  a 
single-frequency  etalon  was  used  instead  of  a 
HeNe  laser.  Its  power  is  variable  up  to  1  W  at  the 
5 14.5-nm  line.  The  detector  amplifier  was  changed 
so  that  the  detector  fed  into  two  separate  ampli¬ 
fiers,  a  low-frequency  amplifier  for  displaying  the 
fringe  position  and  a  wide-bandwidth  amplifier 
(100  kHz  to  50  MHz)  for  displaying  the  acoustic 
pulses.  The  sensitivities  are  1.  6  X  104  V/W  for 
the  wideband  amplifier  and  1.4  X  104  V/W  for 
the  low-frequency  amplifier.  The  main  advantage 
of  this  system  is  that  the  wideband  amplifier  can 
swing  both  positive  and  negative  without  saturating. 
The  experiments  involving  the  side-by-side  gener¬ 
ating  and  detecting  system  were  conducted  to 
characterize  the  signal  generated  on  the  face  of  the 
sample.  The  same  2.5-cm  thick,  5-cm  diameter, 
polished-steel  sample  discussed  in  section  4.2  was 
used.  The  signal  shown  in  figure  16  was  typical  of 
those  obtained.  The  separation  between  the  gener¬ 
ating  and  readout  spots  in  this  case  was  8.6  mm. 
The  ringing  at  the  beginning  of  the  trace  is  electrical 
pickup.  The  first  pulse,  which  is  somewhat  square, 
is  due  to  a  surface  wave  propagating  from  the 
generating  spot  to  the  readout  spot.  The  narrow 
pulse  at  the  center  of  the  picture  corresponds  to  the 
reflection  of  the  bulk  wave  from  the  back  face  of  the 
sample.  The  next  large  peak  is  due  to  reflections  of 
the  surface  waves  from  the  sides  of  the  sample.  A 
flaw  in  the  sample  would  be  visible  somewhere 
before  the  first  narrow  spike  at  the  center  of  the 
picture.  Its  exact  location  on  the  trace  would 
depend,  of  course,  on  the  precise  location  of  the 
flaw,  but  the  reflection  from  the  flaw  would  show 
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Figure  15.  Experimental  setup  for  reflected  acoustic  wave  and  flaw-detection  experiments. 
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Figure  16.  Typical  acoustic  signal  observed  when 
generating  and  readout  spots  are  side  by 
side  on  one  face  of  sample. 


up  as  a  narrow  pulse  before  the  pulse  reflected 
from  the  back  face  of  the  sample.  The  surface  wave 
shown  at  the  beginning  of  the  trace  could  pose  a 
problem  if  the  flaw  in  the  sample  were  positioned 
so  that  its  reflection  would  show  up  at  the  same 
position  as  the  surface  wave. 

Experiments  aimed  at  detecting  a  flaw  in¬ 
side  a  sample  were  carried  out  using  the  steel  sample 
described  above.  To  simulate  a  flaw,  two  2-mm 
diameter  holes  were  drilled  in  the  side  of  the 
sample,  2.0  cm  deep  and  1.65  cm  from  the  front 
face.  The  distance  between  the  holes  was  1  mm. 
The  orientation  of  the  generating  spot  and  the 
readout  spot  with  respect  to  the  flaw  is  shown  in 
figure  17.  From  the  position  of  the  holes  and  the 
velocity  of  the  bulk  wave  in  steel  (about  5.9  X 
105  cm/s),  the  reflection  from  the  holes  should  be 
visible  about  2.7  jis  before  the  reflection  from  the 
back  face  of  the  sample.  Figure  18(a)  shows  the 
signal  received  by  the  readout  system  with  this 
configuration  and  a  generating  laser  power  density 
of  1 40  MW/ cm2.  The  reflected  pulse  from  the  back 
face  is  the  narrow  pulse  near  the  end  of  the  trace, 
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Figure  17. 


Drawing  of  flawed  sample  showing  orientation  of  flaw  and  positions  of  generating  and 
readout  spots. 
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Figure  18.  Oscilloscope  photographs  of  acoustic  waves  observed  in  flawed  sample.  Position  of  flaw  is 
indicated  by  arrows.  Sensitivity  is  increased  by  factor  of  three  in  part  b,  and  5-/xs  delay  is 
added  as  well. 
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and  the  holes  are  indicated  by  the  double  pulse 
visible  about  2.7  /is  earlier.  By  increasing  the 
readout  laser  power  along  with  the  power  output  of 
the  ruby  laser,  the  system  sensitivity  was  increased 
three  times,  as  shown  in  figure  18(b).  Reflections 
from  the  holes  are  easily  visible  in  these  pictures, 
and  the  reflections  can  easily  be  distinguished  from 
each  other.  It  is  evident  that  in  this  case  the  sur¬ 
face  wave  did  not  create  a  problem  for  the  detection 
of  the  interior  flaw. 

Flaw-detection  experiments  directed  at 
observing  flaws  oriented  normal  to  the  surface 
were  conducted  using  a  105-mm  artillery  shell 
casing.  To  accommodate  the  shell  geometry,  the 
readout  system  had  to  be  modified.  The  generating 
laser  beam,  readout  laser  beam,  flaw,  and  shell  po¬ 
sitioning  are  shown  in  figure  19.  Since  the  shell  sur¬ 
face  was  not  polished  and  is  somewhat  diffuse,  a 
focusing  lens  was  needed  to  collect  the  readout 
beam  reflection  from  the  shell.  In  addition,  the 
detection  sensitivity  was  increased  by  the  use  of  an 
avalanche  photodiode  detector.  The  interfero¬ 


meter  was  set  up  with  a  slight  angular  misalign¬ 
ment,  thereby  creating  a  line  or  bar  interference 
pattern  at  the  detector.  A  variable  slit  and  lens  were 
used  to  aperture  this  pattern  and  to  focus  it  onto  the 
detector.  With  these  changes,  the  sensitivity  was 
essentially  the  same  as  when  polished  surfaces  were 
used.  However,  the  mechanical  setup  was  not  as 
stable  as  previously.  The  components  of  the  inter¬ 
ferometer,  including  the  sample,  could  not  all  be 
mounted  securely  to  one  plate  as  with  the  small 
steel  sample  experiments,  and  vibrations  in  the 
system  caused  fringe  shifts .  Also,  due  to  constraints 
in  mounting  the  components ,  the  path  length  in  the 
two  legs  of  the  interferometer  could  not  be  made 
nearly  the  same,  thereby  creating  fringe  motion 
due  to  laser  frequency  variation.  These  problems 
could  be  corrected  in  the  following  manner.  First,  the 
laser  stability  in  both  amplitude  and  frequency  could 
be  improved  by  isolating  the  laser  from  reflections 
from  the  interferometer.  Secondly,  a  servo¬ 
mechanism  could  slave  the  mirror  position  to  the 
fringe  position,  thereby  reducing  the  variation  in 
the  fringe  pattern  caused  by  mechanical  vibra- 
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tions.  These  changes  would  increase  the  readout 
system  sensitivity  and  reproducibility  quite 
markedly.  However,  due  to  time  constraints,  it  was 
decided  to  use  the  system  with  the  former  defects, 
even  though  the  combination  of  the  two  made  it 
very  difficult  to  get  reproducible  data,  as  the 
position  of  the  fringe  pattern  prior  to  the  laser  pulse 
was  random.  The  effect  of  this  can  be  readily  seen 
from  figure  4. 

Flaws  that  originate  in  the  1 05-mm  shells 
during  manufacturing  usually  run  along  the  length 
of  the  shell,  penetrating  into  the  shell  normal  to  the 
shell  surface.  They  occur  most  frequently  on  the 
inside  or  outside  surfaces  of  the  shell  casing  wall. 
Flaws  with  this  orientation  are  inherently  more 
difficult  to  detect  with  this  system  than  internally 
occurring  flaws.  To  simulate  the  geometry  of 
manufacture-generated  flaws,  several  slots  with 
different  dimensions  were  cut  along  the  length  and 
on  the  outer  surface  of  the  shell.  The  slot  sizes  were 
(a)  4  mm  deep,  0.5  mm  wide,  and  10  mm  long;  (b) 


2  mm  deep,  0.5  mm  wide,  and  10  mm  long;  and(c) 
0.5  mm  deep,  0.25  mm  wide,  and  6  mm  long. 
Flaws  (a)  and  (b)  were  detectable,  but  flaw  (c)  was 
not,  probably  due  to  the  instability  in  the  interfero¬ 
meter  discussed  previously.  To  detect  the  flaws, 
the  generating  spot  was  positioned  1 5  mm  from  the 
readout  spot,  with  the  flaw  being  between  the  two, 
but  as  close  to  the  readout  spot  as  possible.  The 
generating  spot  area  was  0.14  cm2  and  the  power 
density  at  the  shell  was  170  MW/cm2.  With  this 
configuration  the  traces  in  figures  20(a)  and  (b) 
were  obtained  for  the  2-mm  deep  flaw.  Figures 
20(c)  and  (d)  show  similar  signals  obtained  in 
unflawed  regions  for  the  same  scale  settings. 
Signals  in  various  unflawed  regions  showed  similar 
pulse  shapes  as  those  shown  in  figures  20(c)  and 
(d).  The  signals  obtained  with  the  readout  spot  in 
an  unflawed  region  are  made  up  of  a  surface  wave 
that  propagates  directly  from  the  generating  spot  to 
the  readout  spot  combined  with  bulk  wave  reflec¬ 
tions  riding  on  top  of  it.  The  bulk  waves  produce 
the  narrow  spikes  that  appear  on  the  slower  surface 
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Figure  20.  Oscilloscope  photographs  of  signal  observed  for  2-mm  deep,  10- mm  long  by  0.5-mm  wide  flaw 
as  compared  to  an  unflawed  position.  The  vertical  scales  are  all  50  mV/div. 
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wave.  The  traces  obtained  with  the  readout  spot 
near  the  flaw  are  clearly  different  than  those  from 
an  unflawed  region.  The  beginning  part  of  the 
signal  is  affected  quite  noticeably  by  the  flaw,  while 
the  later  part  of  the  pulse  is  also  different,  but  not  so 
markedly.  This  size  flaw  is  easily  detectable  by 
comparing  the  initial  part  of  the  pulses,  and  smaller 
flaws  could  be  detected  once  the  various  acoustic 
waves  in  this  geometry  were  identified  and  the 
interferometer  stabilized.  The  signals  observed 
with  flaw  (c)  were  not  significantly  different  than 
those  observed  in  an  unflawed  region.  Since  flaw 
(b)  was  easily  detectable,  the  detection  limit  for 
this  setup  would  be  smaller  than  flaw  (b)  although 
larger  than  flaw  (c).  However,  the  specific  limit  is 
not  known  since  we  did  not  observe  the  signals  for  a 
large  number  of  different  sizes  of  flaws  near  the 
detection  limit. 


5.  CONCLUSIONS 

We  have  demonstrated  that  a  high-power  laser 
pulse  can  generate  a  longitudinal  acoustic  wave  in 
a  metal  sample  suitable  for  use  in  flaw  detection. 
For  a  laser  power  density  greater  than 
50  MW/cm2  we  have  detected  an  acoustic  pulse 
of  the  same  temporal  shape  as  the  laser  pulse  using 
a  BaTi03  transducer.  The  acoustic  pulse  width 
duplicates  the  laser  pulse  width  and  the  amplitude 
makes  the  acoustic  pulse  easily  detectable. 

We  have  also  shown  that  under  a  somewhat 
idealized  situation,  the  sample  under  test  being 
ground  and  polished,  a  laser  interferometer  receiver 
system  can  easily  detect  the  generated  acoustic 
wave.  Signal-to-noise  ratios  for  the  detection  of  the 
transmitted  acoustic  wave  have  approached  50. 
With  an  unpolished  sample  we  were  able  to  get  an 
equivalent  sensitivity  by  modifying  the  readout 
system. 

We  have  also  demonstrated  that  a  side-by-side 
flaw-detection  system,  as  proposed,  can  detect 
flaws  in  solids,  and  is  a  distinct  possibility  for  a 


practical  flaw-detection  system.  We  have  shown 
that  with  a  side-by-side  system  the  acoustic  wave 
reflected  from  the  rear  face  of  the  sample  is  easily 
detectable.  A  simulated  flaw  in  the  sample,  which 
consisted  of  two  2-mm  diameter  holes  drilled  side 
by  side  in  the  sample,  was  easily  detected,  and  the 
holes  were  easily  distinguishable  from  each  other. 
The  practicality  of  the  system  for  production-line 
use  must  still  be  demonstrated. 

Probably  the  most  significant  achievement  was 
the  detection  of  a  longitudinal  flaw  in  a  105-mm 
shell  casing.  We  were  able  to  detect  a  flaw,  2  mm 
deep  by  10  mm  long,  cut  into  the  shell.  This  is  less 
than  one  half  the  size  of  a  critical  flaw  in  this  type  of 
shell  casing.  A  critical  flaw  is  a  flaw  that  has  a  50- 
percent  chance  of  causing  the  shell  to  fail  during 
firing. 

Improvements  in  the  readout  system  would 
allow  the  detection  of  even  smaller  flaws  that  would 
give  a  larger  margin  of  safety.  With  careful  redesign 
of  the  readout  system,  the  stability  problems 
discussed  in  section  4.3  can  be  solved  and  overall 
detection  sensitivity  can  be  improved.  This  redesign 
should  include  optically  isolating  the  laser  from  re¬ 
flections  by  the  interferometer;  mounting  the  fixed 
mirror,  beam  splitter,  and  focusing  lens  on  one 
base;  designing  a  diffraction-limited  focusing  lens; 
using  a  variable  density  optical  filter  in  the  refer¬ 
ence  leg  to  equalize  the  power  in  the  reference  and 
test  legs;  improving  the  detector  electronics;  and, 
finally,  using  a  piezoelectric  drive  to  slave  the 
fixed  mirror  to  maintain  constant  and  optimum 
phase.  The  last  improvement  would  have  the  effect 
of  stabilizing  the  fringe  pattern  for  low-frequency 
drift,  while  not  affecting  the  observation  of  the 
(much  higher  frequency)  acoustic  signals. 

Since  a  critical  flaw  is  quite  large  for  a  1 05-mm 
shell,  being  5  mm  deep  and  about  50  mm  long,  the 
minimal  amount  of  surface  damage  done  to  the 
shell  during  the  acoustic  wave  generation  process 
would  not  create  a  problem.  For  other  applications 
the  amount  of  surface  damage  allowed  would  have 
to  be  determined. 
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